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Abstract We summarize the method of mass modelling of galaxy clusters based on reproduc- 
ing the dispersion and kurtosis of the projected velocity distribution of galaxies. The models 
are parametrized within the framework of the NFW density profile, characterized by the virial 
mass and concentration, together with the constant anisotropy of galaxy orbits. The use of 
velocity dispersion alone does not allow to constrain all the three parameters from kinematic 
data due to the mass-anisotropy degeneracy. The degeneracy is broken by introducing the 
fourth velocity moment, the kurtosis. We tested the method based on fitting both moments on 
mock data sets drawn from simulated dark matter haloes and showed it to reproduce reliably 
the properties of the haloes. The method has been applied to estimate the mass, concentration 
and anisotropy of more than 20 clusters which allowed us to confirm, for the first time using 
kinematic data, the mass-concentration relation found in N-body simulations. 



1 The method of velocity moments 

The traditional approach to studying the kinematic data for galaxy clusters 
involves the modelling of the velocity dispersion profile. The procedure consists 
of collecting enough (preferably more than a hundred) redshifts and positions 
of the galaxies, binning the data and calculating the dispersion values as a 
function of distance from the cluster centre. In the parametric approach one 
then assumes that the mass distribution is given by some model, e.g. the NFW 
(Navarro, Frenk & White 1997) profile characterized by the virial mass M v and 
concentration c, and makes some assumption about the anisotropy of galaxy 
orbits (3 (we assume /3=const with radius). The modelling then involves solving 
the lowest-order Jeans equation for spherical systems (Binney & Tremaine 1987), 
projecting the solutions along the line of sight and adjusting the free parameters. 

An example of such modelling for the case of Abell 576 is shown in Figure 1. 
As we can see, for any anisotropy parameter (3 a value of concentration c of 
the density profile can be found which provides a good fit, which signifies the 
degeneracy between the parameters. Both parameters cannot be constrained 
from the dispersion alone. One can break this degeneracy by considering the 
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Figure 1: The velocity dispersion profile of the cluster Abell 576. The data 
points show the dispersions calculated with 30 galaxies per bin with standard 
sampling errors. The three lines show three examples of projected solutions of 
the Jeans equation with indicated parameters. 



fourth-order velocity moment, the kurtosis, which depends mainly on anisotropy 
(see Sanchis et al. 2004). The improved method thus involves solving two Jeans 
equations for the radial moments (see Lokas 2002) 
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where v is the distribution of the tracer (which we assume to follow the mass), 
I(R) is its 2D projection and <& is the gravitational potential. It is convenient to 
express the fourth moment in the form of the line-of-sight or projected kurtosis 
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Figure 2: Mock kinematic data sets generated from the simulated dark matter 
haloes. 

2 Tests on N-body simulations 

The method has been tested on mock data samples generated from simulated 
dark matter haloes (for details see Sanchis et al. 2004 and Lokas et al. 2006). We 
have used a set of ten dark matter haloes formed in a standard ACDM cosmology 
described by Wojtak et al. (2005). The mock data sets were constructed by 
placing an imaginary observer at a distance of 100 Mpc from a given halo and 
projecting all particle velocities along the observer's line of sight (for simplicity: 
along one of the axes of the simulation box) and their positions on the plane 
of the sky. We then randomly selected 300 particles from the projected phase 
space of each halo limiting the particle velocities to the range ±4000 km s _1 
with respect to the halo mean velocity. 

Three examples of such mock data sets are shown in Figure 2 for one of the 
haloes (halo 2 in Wojtak et al. 2005 and Lokas et al. 2006) with the line of 
sight along the x, y and z axis of the simulation box. Open symbols mark the 
particles unbound to the cluster but present in the samples due to projection 
effects. Such contamination has to be removed before the proper modelling 
at least in the case of the strongest outliers which would otherwise artificially 
inflate the values of velocity moments and bias the estimated parameters of the 
cluster. This can be done by a number of methods (for a review see Wojtak et 
al. 2007). Once the samples are cleaned of interlopers we calculate the profiles 
of velocity moments and fit them with solutions of the Jeans equations (3)- (4) 
estimating the best- fitting parameters M v , c and 0. The results of such tests 
are presented in Figure 3 for all ten haloes observed along the y axis of the 
simulation box. The recovered parameter values in most cases agree within 
errors with the true values measured from the full 3D information about the 
haloes. 
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Figure 3: Comparison of the fitted parameters (dots with error bars) with the 
true parameter values measured from the simulation (crosses) . The virial masses 
are in units of 10 14 M Q . 



3 Application to nearby clusters 

The method of velocity moments as described above has been applied up till now 
to more than 20 nearby (z < 0.1) rich galaxy clusters starting with the Coma 
cluster (Lokas & Mamon 2003), and following-up with other nearby relaxed 
clusters (Lokas et al. 2006; Wojtak & Lokas 2007). An example of the results 
for the cluster Abell 1060 is presented in Figure 4. The left panel shows the 
kinematic sample of 330 member galaxies used for the modelling with two open 
symbols indicating the rejected outliers. In the middle panel wc plot the velocity 
dispersion profile, which as expected for systems with NFW-like density profile 
and isotropic orbits (Lokas & Mamon 2001), decreases slowly with distance. 
The right panel plots the kurtosis-like variable k = (log/*) 1 / 10 which has the 
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Figure 4: The kinematic sample used for the modelling of A1060 (left), the 
velocity dispersion profile (middle) and the kurtosis profile (right). The lines 
show the best-fitting solutions of the Jeans equations. 
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Figure 5: The mass-concentration relation from kinematic data for 22 clusters. 
The solid line is the fit to the data points, the dashed line is the prediction from 
the simulations of Bullock et al. (2001). 



advantage of a Gaussian sampling distribution (see Lokas & Mamon 2003). In 
both panels the solid lines show the best-fitting solutions of the Jeans equations 
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where the errors are \a uncertainties following from the sampling errors of 
velocity moments. 

Similar modelling for other clusters has been performed giving us the pos- 
sibility to verify the mass-concentration relation postulated e.g. by NFW and 
Bullock et al. (2001). The relation obtained by combining data for 22 clus- 
ters is shown in Figure 5. The fit to the data points of the relation logc = 
alog[M„/(10 14 M©)] +b gives a = -0.17 ±0.65 and b = 1.14 ±0.62 (la error 
bars). This fit is plotted as a solid line in Figure 5. We can see that the slope 
of the relation agrees well with the prediction from N-body simulations of dark 
matter haloes by Bullock et al. (2001) shown with the dashed line, although 
the concentrations from the data are typically shifted towards higher values for 
a given mass. 
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